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We apply the renormalization group (RG) method to examine the observable scaling properties 
in Newtonian cosmology. The original scaling properties of the equations of motion in our model 
are modified for averaged observables on constant time slices. In the RG flow diagram, we find 
three robust fixed points: Einstein-de Sitter, Milne and Quiescent fixed points. Their stability (or 
instability) property does not change under the effect of fluctuations. Inspired by the inflationary 
scenario in the early Universe, we set the Einstein-de Sitter fixed point with small fluctuations, 
as the boundary condition at the horizon scale. Solving the RG equations under this boundary 
condition toward the smaller scales, we find generic behavior of observables such that the density 
parameter Q decreases, while the Hubble parameter H increases for smaller averaging volume. The 
quantitative scaling properties are analyzed by calculating the characteristic exponents around each 
fixed point. Finally we argue the possible fractal structure of the Universe beyond the horizon scale. 
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I. INTRODUCTION 



Several observational evidences imply scaling properties in the present Universe. For example, (a) G. de Vaucouleurs 
compiled data of density-size relation for galaxies and clusters of galaxies, and pointed out the scaling relation of 
\ogp — — 1.71ogL, where L is the linear scale of each astrophysical object and p is the mass density averaged at this 
scale 0. This relation asserts the systematic decrease of the density at larger scale, (b) It is widely known that the 
observed two-point correlation function ^(r) for galaxies or for clusters has the scahng property as ^(r) cx r~^-* 
(c) Pietronero and his collaborators discuss a fractal structure of the Universe and claim that the fractal dimension is 
about two up to the scale of 1000 Mpc (d) The observations of mass-luminosity ratio for the various astronomical 
objects linearly increases toward the larger scale up to the scale about lOMpc |4|. 

Motivated by the above observational evidences, we explore systematic underlying physics which may govern the 
scaling properties of observables in the present Universe. Usually the above scaling properties are understood as 
a result of an interplay of the scale invariant initial condition (Harrison-Zel'dovich spectrum) and the gravitational 
instability ||^. In this usual approach, a global uniform background is assumed. In our approach, we do not assume 
such uniform background a priori but consider general inhomogeneous distributions. 

Let us consider the scaling properties carefully. We start our analysis with the equations of motion for fluid 
in the Newtonian Universe neglecting the pressure term for simplicity. This set of equations admits a naive scale 
invariance; a scaling of space-time coordinates as well as an appropriate scaling of observables leaves the set of 
equations unchanged. We define this naive scaling transformation as the operation-S. Such a naive scaling property is, 
however, not necessarily reflected in the actual observations in a direct form. This is because in general the physically 
measurable quantities are strongly affected by the actual method of observations. Therefore the above naive scaling 
should be modified when we apply it to the actual observations. There are at least two important factors for the 
modification: (1) An observable region is limited by causality. (2) We may merely observe averaged quantities in a 
certain region in space. 

We first consider the factor (1). Actually, we can only observe a light signal from a galaxy located on our past 
light-cone. In Newtonian cosmology, the observable quantities are located on a hyper-slice of constant time. Because 
in general the above naive scale transformation also requires a change of the time variable, we have to adjust it so 
that we can compare the observable quantities with different scales on the same time slice. This adjustment can be 
implemented by the time evolution of the observables by using the equations of motion for them. We define this 
adjustment as the operation-T. In this context, we have to reconsider exactly what we are measuring. In this paper 
we consider scaling behavior of the cosmological parameters. For example, the density parameter 17 is a dimensionless 
quantity and therefore it is invariant under the naive scale transformation. However in many cases, we can directly 
measure the matter density p and separately the cosmic expansion parameter H The combination of these 
observations induces the value of the density parameter ^ which depends on the scale. In this case, the scale 
transformation and the necessary time adjustment are non-trivial. Thus we have to carefully distinguish the direct 
observables and the secondary induced observables. 

Next we consider the factor (2). What we observe is necessarily an averaged quantity on some region of space- 
time because there are always spatial fluctuations and the resolution of our observation is limited. Therefore the 
observables are most generally dependent on the scale of averaging. We define this averaging as the operation-A. 
However, it would be very difficult to consider the invariant averaging procedure in the full general relativity and to 
make connection with observables 0. Therefore mainly from this reason we would better restrict our considerations 
to the Newtonian cosmology, in which an averaging procedure is properly defined and is explicitly calculated 

After the above adjustments (1) the operation-T and (2) the operation-A on the original naive scaling properties of 
the underlying equations of motion, we would obtain the relevant information for actually observable quantities. Thus 
in general, the difference of two observables at different scales would be related with each other by the naive scaling 
(operation-S) compensating time evolution (operation-T) and the inevitable averaging procedure (operation-A). 

In order to unify the above operations, S, T, and A, we need to introduce one more ingredient in our considerations. 
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The most sophisticated method so far to deal with such multiple operations would be the renormalization group (RG) 
method §. This is a very general method to obtain the observable response of the system against the scale change and 
is widely used in various fields in physics including quantum field theory, statistical mechanics, fluid mechanics, etc. 
For example in the quantum field theory, a naive scaling in the classical theory should be modified by the quantum 
fluctuations. In the statistical physics, the scaling property is used with the averaging operation to yield critical 
exponents. 

In this paper, we try to apply this RG method to Newtonian Cosmology. In this model, the scaling invariance holds 
all the time. Thanks to the averaging procedure, we can define classical fluctuations of the observables. Then, naive 
scaling behaviors of the observables are modified by renormalization of the classical fluctuations. This is conceptually 
the same as RG methods in other fields of physics. 

The RG equations we obtain give the effective scale change of averaged observables, and the effect of spatial 
fluctuations on observables. The flow diagram generated by the RG equations in the parameter space of observables 
represents how the observables actually change under the scale transformation. Each flow line represents a single 
physical system, in which the averaged variable at different scale corresponds to each point on the flow. Thus the 
RG flow gives a whole set of possible cosmological models. The fixed points, i.e. the stagnation points in the flow 
diagram, of the RG equations represent the possible "uniform" structure of observables, that is the averaged variables 
do not depend on the scale, and characterize diversities of the system. 

There have been some applications of the RG method in cosmology so far ||l^. The RG method is an effective 
tool to examine the asymptotic behavior of a dynamical system. In fact, Koike, Hara and Adachi succeeded in 
explaining the critical behavior of a self-similar solution during the gravitational collapse of a radiation fluid using 
the RG method [0. Subsequent series of studies are in this line of study |l2j. In their cases, the aspects of (S) and 
(T) are reflected in their method. Thus their RG method does not necessarily contain an averaging procedure (A). 
On the other hand, Carfora and Marzuoli and Carfora and Piotrkowska |Q explicitly accounted for the spatial 
fluctuation of physical quantities in general relativity and derived the RG equation reflecting the aspects (S) and (A). 
Since they concentrated on the scaling in the spatial direction, the aspect (T) has not been manifestly reflected in 
their approach. Here in this paper, to examine the scale dependence of the observables on a constant hyper-slice, 
we apply the RG method from the full aspects of (S),(T), and (A) altogether, assuming some scaling property in 
averaging observables. 

This paper is constructed as follows: In section II, we derive the RG equation for the averaged Newtonian fluid. 
The section III is the main part of our analysis. Here we examine the RG flow in the parameter space and show 
that there appear three robust flxed points including the Einstcin-dc Sitter universe. We examine the effects of shear, 
tidal force, and fluctuations, and derive the scale dependence of the physical variables such as 17 and H around the 
Einstein-de Sitter universe. In section IV, we summarize our results. We also present the full RG equations including 
fluctuations up to the second order cumulants in the Appendix. 



II. DERIVATION OF RENORMALIZATION GROUP EQUATIONS 



In Newtonian cosmology, the following equations describe the evolution of a self-gravitating fluid in the Eulerian 
coordinate system: 

I — <-)• 

— = -{v ■ V)v + g , 
V • g = -AirGp , (1) 

where p, v, and g (= —Vcf) respectively denote the mass density, the velocity of the fluid, and the gravitational 
acceleration. Now we transform these equations into the Lagrangian coordinate system. In order to eliminate the 
gravitational acceleration g, we introduce the spatial derivative of the velocity flcld and decompose it into a trace 
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(2) 



part (expansion 6 = V ■ v), a symmetric traceless part (shear cr*^) and an anti-symmetric part (rotation w*^ ) ( 
Then the above equations become 

dp . 

'E = J_e^ + 2{u^-a^)-A^Gp, (3) 

^ = -^^^^ ~ - ^l^-"- "')^^ - ^'r (4) 

<iw\- 2 ■ ■ I. J, 

^ = -/-^-^>'- (5) 
where the Lagrange fuU time-derivative 

d d „ 

is used on the left hand side. In these equations, a and uj denote the magnitudes of shear and rotation, respectively, 
defined by 

^' = ^^>/ and u;2 = iwVw/'. (7) 

The tensor ^ — V*Vj0— (5*j-V^0/3 denotes the tidal force which comes from the spatial difference of the gravitational 
acceleration. In pure Newtonian dynamics, this is a slaving variable and there is no equation of motion for it. 
Fortunately, several excellent approximation schemes have been developed to deal with this highly non-local quantity. 
In this paper, we mainly follow the local tidal approximation (LTA) jist JlGf , in which the equation of motion for the 
tensor E"^ ^ is given by 

^-0 E'j ~ 4ttG p a'j . (8) 

Hereafter we neglect the rotation for simplicity (w*^- = ). Actually the initial assumption of w'^ = is sufficient 
to guarantee this condition all the time as is seen from the above equation. Thus irrotational motions form a closed 
subclass of the full dynamics. Then, we can simultaneously diagonalize both the matrices a^j and E^j. We work in 
the frame where they are diagonalized as tXi and Ei (i = 1,2,3), respectively. Since cr'^- and E^^ are traceless, only 
the two components of each are independent. Then it may be sometimes more convenient to introduce the following 
variables a± and E±: 

(cTi, 0-2, 0-3) = (cr+ + \/3cr_, (7+ - V3cr_, -2cr+), 
{Ei,E2,E3) ^ {E+ + V3E-,E+ - V3E-,-2E+). (9) 

At a glance, we observe a scaling property in the above set of equations of motion (§) - (^), (^. They are invariant 
under the following scaling transformation: 



t - 


^t' = 




X - 


^x' = 


e^^x 


9- 


^9' = 




cr± - 




- e a± 


E± - 




= e-"^^E± 


P - 







(10) 



where s is a constant free parameter and At = t' — t measures a change of scale. 
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Now, based on these scaling properties, we derive renormalization group equations for spatially averaged observables 
on a constant time slice. To analyze the structure of the Universe, we adopt a scaling solution of our basic equations. 
Such a scaling solution could be realized as a result of a self-gravitating dynamical system with a scale invariant 
density fluctuation B, although here wc do not specify any models to give such a scaling solution. If, however, we 



consider a scaling solution as (10) as it is, we find either a homogeneous Universe or inhomogeneous Universe with a 
"center", which is not consistent with the cosmological principle. We are interested in the inhomogeneous Universe 
model consistent with the cosmological principle such as a fractal universe. If we have a fractal universe, we expect 
that some averaged values around an observer may show a scaling property. Such a scaling property will be always 
observed anywhere because of a fractal structure of the Universe, which is consistent with the cosmological principle. 
This is why we consider spatially averaged observables here. We assume a scaling property as ( p^ ) not only for 
averaged values of the variables in the above equations but also for those of the fluctuations such as a second order 
cumulant. 

We define a spatial average (...) of an observable /(t, x) at time t as 

= T777T / d'^f{t,x), (11) 



Vit) 



V(t) 



where Vit) is a spatial domain at time t and V{t) is its volume defined by V{t) — jx>(t) '^^^ observable /(t, x) 
can be any function of 6, p, a±, E±. We would like to derive the expression for the infinitesimal scale change of the 
averaged observable on the same time slice, i.e. 

^{f)v(t) _ {f)v>(t) - {f)v(t) , . 

Ai ~ Al ' ^ ' 

where the parameter £ measures the true scale change on this time slice, which is defined by 



(v'it)\ 



1/3 



»^=»*-'=(vWj ■ '''' 

and is different from the naive scaling parameter t. "D' [t) and V {t) denote the domain and volume transformed back 
to the original time t-slice by time evolution from 'D'{t') and V'{t'), which are obtained by a scale transformation 
(|l^). The above equation (|l|) can be decomposed into two parts: 

(i) a change As{f) associated with a scale transformation from I?(t) to 'D'{t'), 

(ii) a change Axif) associated with a time evolution from 'D'{t') to ^'{t). 
That is, 

= ^ [{{f)v'{t') - {f)v(t)) + {{f)v'(t) - (/)p'(t'))] 

^ ^^{AsU) + Arif)) . (14) 

As for part (i), we use the scaling law of the observable {f)T>{t), which is 

{f)vit)^{f)v'it')^e'^^^^f)j,^,y (15) 

The second part (ii) is given by the Taylor expansion of f{t',x') and 'D'{t') with respect to time. 

1 



f{t,x') + f^{t,x')At 



0{{AtY) 



{f)v'it) + M l-{9)j,,^t){f)v'it) + mv'it) + l^f^^^ ^ j + 0{{Atf), (16) 



where At = t' — t. 
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The true scale change A£ on a constant time shce is given as fohows. In the scale transformation (i), the volume 
change is V {t')/V{t) = e^'^'^ = 1 + 3Z\t, but it includes the expansion effect of the Universe because time is also 
transformed. To obtain the true scale change on the same time slice, we have to transform it back to the original 
time slice by the time evolution (ii), which gives V' [t') /V' [t) = 1 + {9)At. Then we have 



\vit)J 



= 1 
= 1 



Ai. 



Thus we obtain 



The time interval At = t' 



AI = At- ^At. 
3 



t is given by the scale transformation (10) as 
At 

= St. 



At 

Hereafter we identify this fixed time t as the present cosmic time to, when we observe cosmological quantities. 
Thus, combining Eqs.(14) ^ (p^, we obtain the differential equation for the averaged observable {f)v(ty- 



(17) 



(18) 



(19) 



d{f) 



V{t) 



di 



S\-{f)v(t) + {0)vit){f)vit)~{Of)v(t) 




where 



S 



(20) 



(21) 



l-sWp(t)^o/3' 

The parameter s-dependence appears solely through this factor S. 

The equation ( po| ) is not yet sufficient to obtain a closed set of RG equations because of the non-linearity of the 
basic equations. Differential equations for the n-th order averaged quantities are not closed up to the n-th order. 
They necessarily contain {n + l)-th order averaged quantities. This is the famous BBGKY hierarchy in statistical 
physics. The ordinary method to obtain the closed set of equations is to truncate this hierarchy at some order. In our 
case, we include the effect of fluctuations at the lowest non-trivial level, i.e. the second order cumulants. Then we will 
neglect all intrinsic fluctuations (cumulants) higher than the second order, finding the following truncation formula: 



ifgh) ^ {f){9h) + {g){fh) + {h){fg) - 2{f){g){h) 



(22) 



To write down the basic equations, we shall use the second order cumulant {fg)c instead of the averaged quadratic 
quantity (fg), where {fg)c = {fg) — {f){9)- Applying this reduction, we obtain a closed set of 27 differential equations 
for averaged variables and their second order cumulants. We will present the complete expression in the Appendix. 
Here we only show a reduced set of equations to the first order (the averaged variables) by setting {fg)c ~ 0: 

d{p) 



de 

m 

di 
de 
de 

d{E+) 

de 

d{E-) 



= s 



(0) 



3 



(P>, 

1 I {9)+A7rG{p)+6i{a+f + {a^f) 



1 



-(0)-l)(a_)+2(a+)(a_) + (£;_) 



(23) 
(24) 
(25) 
(26) 
(27) 
(28) 
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These also form a closed set of equations at this order. We call these the first order RG equations, while those in the 
Appendix the second order RG equations. 

Here and in Appendix, we have set to — 1 and dropped the sufhx of integration domain for averaged values just 
for simplicity. Hence, in this units, the expansion 6 and the density p for Einstein-de Sitter universe turn out to be 
0Eds = 2 and pEds = 1/QttG 



III. ANALYSIS OF RG EQUATIONS : FIXED POINTS AND RG FLOW 

We now examine the RG equations obtained in §. 2. First, we analyze the first order RG equations ( p3| ) ^ ([2^). 
Although this may not be realistic because the structure of the Universe could be highly inhomogeneous in our model, 
it shows a naive idea to understand the RG flow and give three important fixed points (including the Einstein-de 
Sitter spacetime), which will be also most relevant even in the case with fluctuations. We then examine the stability 



of these fixed points Next, introducing fluctuations, we study the second order RG equations ( |AlD ~ ( |A27| ) 

and analyze the stability of those fixed points in a much wider parameter space. Finally in this section, we set the 
Einstein-de Sitter fixed point with small amount of fluctuations, as the boundary condition at the horizon scale in 
our present Universe and examine a scale dependence of averaged variables toward smaller scale. 



A. Fixed points of the RG equations and Stability 



First we analyze the structure of our RG equations. From Eqs. (Al) ^ (A27), we find that those are written in a 
vector form as 



d{F) 

de 

d{FF% 

de 



^S{A[{F)]+C-{FF),) 
= S B[(F)]-(FF),, 



(29) 
(30) 



where (F) and {FF)c are 6- and 21-dimcnsional vectors defined as (F) = {{fi}) and {FF)c = {{fifj)c) with fi's 
(i = 1 ~ 6) denoting p, 9, <t^, (7_, and E_. A[{F)], B[{F)], and C are 6-dimensional vector of quadratic functions 
of (/j), 21x21 matrix of linear functions of (fi), and 6x21 constant matrix, respectively. Setting {FF)c — 0, we have 



d{F) 

de 



SA[{F)], 



(31) 



which is exactly the first order RG equations (|2^) ~ (|2§|). 

To analyze the structure of the RG equations, we first need to know the fixed point ((F), (FF)) = {{F)* , (FF)*), 
which is defined by 



d{F) 



de 



d{FF), 



{F) = {F)*,{FF)^ = {FF)* 



de 



0. 



(F) = (F)%(_FF)e = (-FF); 



As for the first order RG equations, the fixed point (F) — (F)* is given by 

A[{F)*] = 0. 



(32) 



(33) 



The above equations (p9|),(|30|) with ( p3| ) guarantee that the fixed points (F)* in the first order RG equations are 
always those for the second order RG equations with (FF)* = 0, although new additional fixed points appear in the 
second order. 

The stability analysis is very important to know asymptotic behaviors near the fixed points. The stability is 
examined by linearizing the RG equations around the fixed point (F)* , (FF)* as 
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{F) = {F)*+T, 



(34) 
(35) 



where both and {J-J-)c are small enough to make the linear perturbation treatment effective around (i^)*, (FF)*^ 
Then, we find the perturbation equations in the form 




Wll Wl2 

W21 W22 




(36) 



where the perturbation matrix Wl^l = (Wab) is given as 



9{SA[{F)]) 
d{F) 



d{SBl{F)]) 
9{F) 



{F) = {Fy 



Wll W12 
>V2i W22 

\ d(F) 

{F) = {F) 

with S* — S\(^p^^(^py . For the first order RG equations, it turns out to be 



S*C ^ 
S*B[{Fr] ^ 



(37) 



where Wl^l is given as 



,dA[{F)] 



d{F) 



(38) 



(39) 



{F) = {F)' 



because of Eq. (^ For the fixed points in the first order RG equations, we find the following important results. As 
we mentioned, they are also the fixed points in the second order with (FF)* = 0. Then the perturbation matrix Wl^l 
becomes 






(40) 



This form guarantees that the eigenvalues of the first order matrix W^^l are always those of the second order matrix 
yyl^l . As we will show later, the other eigenvalues of the second order matrix are constructed by a sum of any two 
eigenvalues of the first order matrix (see Eqs. (|48|), (^l|)). 

The positive (negative) eigenvalue of the matrix W^'^l (or W^^l) shows the instability (stability) of the fixed point 
{{F)* , (FF)*) (or (F)*) toward the larger scale, i.e. in the large ^-direction. The stability and the instability are 
exchanged if we follow the RG flow in the opposite direction, i.e. toward the smaller scale. 

The eigenvalues fully characterize the scale dependence of all physical variables in the vicinity of each fixed point. 
Although each observable is scaled according to its dimension, the effect of dynamics changes the scaling property 
from the original bare scaling law (|l^), so that the scaling of each observable is determined by the superposition of 
several eigen functions as follows: 



(/.) - (/.)* + 

Amin = -^1 < -^2 < 



N 

771 — 1 



• < Aat = A 



(41) 



where (/i)* is the value of (fi) at a fixed point, e^^m is some small constant, and Am {m = 1 ~ N) are eigenvalues. 
N = n{n + 3)/2 for the second order RG equations, while = n for the first order RG equations, with n being the 
number of observables. Since the maximum (minimum) eigenvalue, Amax(Arnin) becomes the most relevant toward 
larger (smaller) scale, the scaling property of each observable around fixed point is determined by Amax(Arnin)- 

The fixed points themselves are independent of the scaling parameter s, but the matrix W'^l (or W'^') and its 
eigenvalues depend on s. For the sake of explicit argument below, we restrict the range of the parameter s as 
< s < 1 where the topology of the RG flow is not changed |19|. 
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B. Fixed points and the flow of RG equations without fluctuations 



Here we analyze the first order RG equations (23) ~ (g8|) in detail. First, by neglecting both shear and tidal force 
terms, we have simple RG equations ( ^T|) with 

A[{F)] = ^ ' ' \ (42) 

for two averaged variables {F) ~ ((/?), {9)). We easily find three fixed points {k = 1,2,3), which are listed in 

Table 1. 

The stability at each fixed point {F)*f.-^ is examined by linearizing the RG equations around the fixed point (Eq. 
(^)). All eigenvalues of Wl^l for three fixed points are also given in Table 1. For the first fixed point (1) where 
both and vanish, the dynamical term At{F) does not affect the scaling, so that the eigen vector points 

the direction of each observable. On the other hand, for other fixed points, the effect of dynamics changes a scaling 
property from the original bare scaling law ( p^ ) . 

As seen from Table 1, the fixed point (1) is stable toward larger scale. For increasing scale, the Universe approaches 
this fixed point, at which both the energy density and expansion vanish. Therefore we shall call it the fixed point Q as 
the Universe is quiescent. The second fixed point (2) is a saddle point. It represents the "Einstein-de Sitter" Universe 
because {0)*2) — ^EdS and (p)(2) = PEdS- We then call it the fixed point E. The third fixed point (3) is unstable in 
any direction. At first sight, it might look strange since the expansion of the Universe is finite in spite of vanishing 
energy density. This corresponds to the "Milne" Universe in which the expansion rate is apparently finite because of 
a specific choice of the time slice. We call it the fixed point M. This Universe is reduced to the "Minkowski" space by 
adjusting the time slice appropriately. 

In Fig.l, to see more precisely those fixed points and the behaviors of global structure of the RG equations, we 
depict the RG flow on the H-Vl plane, where the Hubble parameter H and the density parameter Q. are defined by 

^ f (43) 

n . (44) 

The RG flow either approaches fixed point Q or escapes to infinity toward larger scale. These two kinds of fiows are 
clearly separated by the divide which passes through two fixed points E and M (Fig.l). The fiow escaping to infinity 
seems unphysical since both H and O blow up to infinity at large scale. 



Next we include shear and tidal force, which leads to Eqs.(23) ~ (|28|). The three fixed points Q, M and E found 
in the case without shear and tidal force survive, and their stabilities also remain unchanged (see Table 2). In this 
sense, the three fixed points Q, M and E are robust. Four new fixed points appear in addition to those three fixed 
points, which are listed in Table 2. Since all new fixed points are saddle points, any RG flow still approaches the fixed 
point Q or escapes to infinity toward larger scale, or it approaches the fixed point M toward smaller scale. 



C. Fixed points and the flow and of RG equations with fluctuations 



Since our model of the Universe may be quite inhomogeneous and could be fractal in a microscopic point of view, an 
inclusion of fluctuations of the physical variables is inevitably important, when we compare our results with those for 
the real Universe. In fact, if we do not include fluctuations, it is mathematically equivalent to a self-similar solution 
on a fixed time slice ||2^, although the physical meaning is quite different. Then, we shall add fiuctuations and see 
how the fluctuations affect the RG flow diagram. We just include second order cumulants to the RG equations and 
truncate this hierarchy there by ignoring all intrinsic fluctuations (cumulants) higher than the second order, although 
we may have another choice of higher-order cumulants [p3| . 
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There are 27 independent variables in this case and it becomes more difficult to find a full list of fixed points. Then, 
first we reduce the number of variables and start our study by the shear-free case. Since independent observables are 
just p and 9 in this case, the total number of variables including second order cumulants is just five: (p), {9), (p^)c, 
((?^)c, {p9)c- We list up all fixed points and show their stabilities in Table 3. The three important fixed points (1)^(3) 
in the first order RG equations still remain and their stabilities are also unchanged. They are quite robust. Since 
all of the new fixed points (8)~(10) are saddle points, global properties of RG fiows toward both larger and smaller 
scales are solely determined by the fixed points Q and M, respectively. The fixed point E remains as a saddle point. 

As we mentioned, the fixed points found in the first order equations are also those in the second order with zero 
cumulants. For the eigen values at those fixed points, we find eigen values in the first order are also those in the 
second order. To see more detail, we explicitly write down the perturbation matrix Wl^l (Eol) as 



\w2i W22 ) \ 47rG l{9)*-l) \^0-|0 

/2((0)*-2) 2{p)* \ /2«;n 2w^2 

S*B[{F)*]^S* { 2(|(0)*-1) %ttG = 2w22 2w2i | (45) 

V AttG " {pY f(0)*-3/ \ W2l Wi2 WII+W22 

Then, the eigenvalue equation for the second order perturbation matrix Wl^l is 

(A^ - trWW • A + detW[il)(A - tr>vW)(A2 - 2 trW^^l • A + 4 detwl^l) = 0. (46) 
Since the eigenvalue equation for the first order matrix Wl^l is 

A^ - trWl^l • A + detWl^l = 0, (47) 
supposing that A^j^'^' and Aj^' are its two solutions, we find the solutions of Eq. (^), Ap' {i = I ^ 5), as 

^[2] ^ ^[1] ^[2] ^ ^[1] 

a[^U2aW, = trWW^AW+A™, A^^l = 2aW. (48) 

Therefore the minimum eigenvalues at fixed points in the first order, in particular the fixed points Q, E and M, are 
given by 

•^min ~ ^'^Lin' (^9) 

if a|j^!„ < (this is the case for the fixed points Q and E), otherwise 

^min — ^mim \^'^) 

(this is the case for the fixed point M). 

Finally, we study a full model with 27 RG equations for 6 independent variables including shear and tidal force and 
their second order cumulants. Though it is very difficult to find all fixed points, all seven fixed points in the first order 
RG equations (listed in Table 2) remain as fixed points, as we mentioned. As for their eigenvalues, we also find the 
same result as that in the case without shear and tidal force, that is, the eigenvalue in the second order RG equations 
{Am \m — 1, ■ ■ ■ , n{n + 3)/2} are constructed by those in the first order {Am |to = 1, • • • , n } as 

A^l = aW for m = l,..-,n 

AI^I = a{^1 + a|,^1 for m = 71 + 1,- ••,n(n + 3)/2 with = I,---, 71 (51) 
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where n is the number of observables. This result is understood from some special relation between the perturbation 



matrices Wl^l and S*B[{F)*] |24|; 



S*BijM[{F)*] = ^ [SikWji + SiiWjk + SjkWii + SjiWik] , (52) 

where fc, I denote the observables p, 0, a±, E±. 

In Table 4, we have only listed three important fixed points Q, E and M with their stabilities. These three fixed 
points still survive and their stabilities remain unchanged even after including fluctuations as well as shear and tidal 
force. 

Let us just mention how our results change if we use another approximation scheme for the equation of motion 
for Eij. We mainly use LTA approximation in this paper. In the case of other local approximations like NMA 
psf , the RG equations become more complicated and it is more difficult to look for all fixed points and eigenvalues. 
However, it is still easy to show that the above three fixed points Q, M and E survive since the RG equations exactly 
reduce to the shear-free case when we set E± = a± = 0. The stabilities of these fixed points can also be checked 
by examining the RG flow around them. We show that the stabilities of these fixed points remain unchanged even 
in the NMA. The minimum eigenvalue, which is the most relevant toward smaller scale, is also invariant irrespective 
of local approximation scheme of Eij, because the eigen-direction for a|^|j^ (or xl^\^^) is orthogonal to both shear and 
tidal directions in 27-dimensional "phase" space of {{F), {FF)c) (or in 6-dimensional "phase" space of (F)). Hence 
the effect of shear or tidal force seems irrelevant to the scaling property of observables around fixed point E or M 
toward smaller scale. 

The three fixed points E, M, and Q are robust and change neither their stability nor the scaling property of 
observables around them toward smaller scale, regardless of shear, tidal force, and then the approximation scheme of 
the tidal force. 



D. The flow of our Universe 

The RG flow and the fixed points we argued in the previous subsections describe all the possible diversity of our 
system. Only a part of them, i.e. a particular integral line in the flow diagram which satisfies certain boundary 
condition, actually describes our real Universe. In order to elucidate such integral flow, we need a plausible boundary 
condition and examine the global/local scaling properties of observables. 

One plausible boundary condition at the present horizon scale and beyond could be the fixed point E ( "Einstein-de 
Sitter" Universe). This boundary condition would be inferred from the infiationary scenario in the early Universe, 
which predicts the flat FRW Universe with tiny fluctuations as our real Universe. 

With this boundary condition, we solve the RG equations toward smaller scale and study the scale change of various 
observables. In solving the RG equations toward the smaller scale, we must remember that (i) the positive (negative) 
eigenvalue of the matrix Wl^l at a fixed point denotes stability (instability), (ii) The direction of the integration is 
opposite to that of flow indicated in Fig. 1. It is clear that the choice of the exact values of point E as a boundary 
condition is trivial and not interesting because it is the flxed point and the flow stays there forever. To find more 
interesting model of the Universe, we need to start from the close vicinity of the fixed point E. 

Here we choose the boundary condition as follows. For the averaged variables, we use the exact value of the fixed 
point E: 

(9) - 0Eds = 2, (p) = pEds = l/67rG, (a±) = (^±) = 0, (53) 

which gives f2 = 1. Since we expect some fluctuations even in an inflationary model, we set small finite fluctuations 
at the fixed E point. We assume that the off-diagonal fluctuations have initially no intrinsic fluctuations ||2^, i.e. 

{f^fJ)c = (z ^ j) (at £ = 0). (54) 
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Before going to show our numerical result, we shall discuss the asymptotic behaviors near three robust fixed points, 
Q, E, and M. We showed that the scaling property of observables around those fixed points toward smaller scale is 

[2] 

determined by the minimum eigenvalue X[-^^^y This is always true for the second order cumulants {fifj)c- The scaling 
property of the averaged variable (fi) may, however, depend on the case, because the eigen vector with the eigenvalue 
'^min perpendicular to the hypcrsurfacc of the cumulants in the "phase" space {{F), {FF)c)- If we have fluctuations 
initially rather than a deviation from the fixed point, i.e. {ff)* 7^ 0, then 

(/»>^(/«>*+e.e^!l^ (55) 

as £ decreases, unless {fi)*c^{{fi)*Y- This is just because how to leave from the fixed point is determined by the 
fluctuations {ff)1, which scaling property is fixed by a|^|j^. On the other hand, if we have an initial deviation from 
the fixed point as {fi) 7^ {fi)* without the initial cumulant ((/f)c — 0), then {fi) — {fi)* behaves first as exp(A|^|j^£). 
Since our boundary condition is the former case, we expect the asymptotic behavior as Eq. (p5|). 

It may be more comprehensive to show our results if we introduce the "dispersions" Aij instead of the cumulants 
{fifj)c- For the diagonal cumulants {ff)c, we define them as 

K = ip')c/{p)', = {e')J{9)^ Al^ = (4)c/W^ = {El)c/{e)\ (56) 

where for shear and tidal force, A^^ and A^;^, which are normalized by {9), have been used rather than the conven- 
tional dispersions because those averaged variables sometimes vanish even if their cumulants are finite, resulting in a 
divergence of the dispersions. The dispersion parameters Ap, Ag, Aa±, Ae± are assumed to be sufficiently small at 
the boundary £ — 0. In what follows, just for simplicity, we also set Ap — Ag = Aq at ^ = 0, which may be plausible 
because p and may be closely connected with each other through the fluid equations and are expected to have the 
same order of fluctuations. 



1. The case without shear and tidal force 

First we consider the case without shear and tidal force. In Fig. 2, we depict the scale dependence of f2, H, Ap, Ag 
in terms of L = e^, which is the ratio of the averaging scale to the horizon scale (or the scale at our boundary). Since 
we have chosen ^ = as the horizon scale (or beyond), our observable Universe is described by negative £. As seen 
from Fig. 2, the integrated flow line always converges to the values at fixed point M, i.e. r^Miinc — and -ffMiino = f • 
This means that the "Milne" Universe is the unique stable fixed point toward smaller scale. Once it comes close 
to fixed point M, the behavior of the RG flow at small scale is determined by the minimum eigenvalue at M, i.e. 

ASn=V(l-^)- 

To see a role of the scaling parameter s, in Fig. 3, we plot the scale dependence of the density parameter fl, the 
expansion H, and the fluctuations Ap, Ag for several values of s with flxed initial fluctuations Aq = 10^^. We flnd 
that the observables show the following universal behavior: 

(i) H increases toward smaller scale from 2/3 to 1. 

(ii) O reduces from 1 to toward smaller scale. There is a plateau near the fixed point E, while a constant slope is 
found near the fixed point M. 

The asymptotic values are given by those at the fixed points E and M. We can also explain a constant slope at the 
smaller scale and a plateau at larger scale in ft diagram. Since {p) ^ 0(exp{[s/(l — s)]£}) and {9) ^ 3 near the flxed 
point M SI at smaller scale behaves as 

l]^0(exp{[s/(l-s)]n), (57) 

where s/(l — s) is the minimum eigenvalue evaluated at the flxed point M. Then the slope of the fi-curve near M 
turns out to be 

dhgn s 
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On the other hand, the plateau appears near the fixed point E, because 



(p)/pEdS ^ 1 + ep exp{-[2s/(3 - 2s)]^}, (0)/0EdS ^ 1 + ee exp{-[2s/(3 - 2s)Y} 



(59) 



and then 



17 - 1 + (ep - 2ee) exp[-(2s/(3 - 2s))£] « 1, 



(60) 



where we have used \y^i^ — — 2s/(3 — 2s) evaluated at E. The turning point Lex of the slope into the plateau is 
estimated easily. Since the deviations from the EdS values are due to the fluctuations Aq, ep,ee ~ 0(Ao). Then, 



we find a large deviation from the EdS values. Thus the minima of the eigen values at two fixed points E and M 
determine the asymptotic structure of our Universe. 

Is there any observational data which may suggest a scale dependence of either f2 or density? In discussion of dark 
matter, when we plot the mass-luminosity ratio in terms of the scale of the objects, we find a similar curve to Fig. 
3(a). If we assume that the mass-luminosity ratio curve is what we are discussing here, we find some constraint of s 
as 0.3 < s <0.7. 



Next, we examine the effect of shear and tidal force. Just for simplicity, we assume that cr_ = E- = 0, which is 
consistent with the RG equations. As we assumed, ((T+) and (£'+) vanish at £ = 0. Then initial free parameters are 
Aq, Ao._.o and KE_fl, which are fiuctuations at £ = 0. Our numerical calculations, shown in Figs. 4, 5, reveal that 
they tend to prevent H and Vl from converging to the fixed point M though "Milne" Universe is still a stable fixed 
point. More precisely, we observe three qualitatively different behaviors of the RG flow of our Universe: 

(i) The RG flow monotonically converges to fixed point M. H and 17 approach to the value determined by fixed 
point M toward smaller scale. 

(ii) Initially the Universe comes near the fixed M point, resulting that H increases and VI decreases toward smaller 
scale. However it turns around before reaching the fixed point M and eventually diverges to infinity in the flow 
diagram as iJ — > — oo and f2 —> oo. 

(iii) The RG flow monotonically diverges to infinity in the fiow diagram as _ff — > — oo and fl oo. 

We can classify the parameter space (Aq, Aa-+fi, Ae+.o) of fiuctuations at the boundary (near the horizon) into three 
regions according to the above asymptotic behaviors of the RG fiow. The class (i) is realized if the fluctuations of shear 
and tidal force are sufficiently smaller than those of the density and the expansion, i.e. Acr+,o/^o ~ ^-E+,o/^o < Rcr- 
Otherwise, i.e. if Ao-^a/Aq ^ Ae^a/Aq > Rcr, the class (iii) will happen. The critical value Rcr is almost 
independent of Aq and approximately Rcr ~ 0.15. The class (ii) is a narrow boundary between them. Hence, 
the class (i) (or (ii)) naturally occurs as long as the shear and tidal force are sufficiently smaller than the density 
fiuctuations at the horizon scale. This situation is strongly suggested by recent observations i.e. the shear is 
strictly restricted around a/H < 10~^ by CMB data. Hence we expect that our Universe belongs to the class (i) and 
its flow line converges into the fixed point M in smaller scale. The s dependence of fl or H is quite similar to the case 
without shear and tidal force (Fig. 6 (a), (b)). 

As for the shear and tidal force, their averaged values increase toward smaller scale, but eventually decrease if the 
flow will converge to the fixed point M (Fig. 6 (c), (d)). The scaling properties of the shear {a+) and the tidal force 



when Aoexp{-[2s/(3- 2s)]4i.} - 0(1), i.e. 



InLcr = 4r ^ [(3 - 2s)/2s] In Ao 



(61) 



2. The case with shear and tidal force 
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are also determined by Ajjjjj at the fixed points E and M. Both shear and tidal force are dumped off at smaller 
scale along the flow converging to the fixed point M. 

The behaviors of fluctuations Ap, Ag, Ao-_|_ and A^;^ are shown in Fig. 6 (e)~(h). We find that the similar behaviors 
to the case without shear and tidal force for Ap and A^i. 

E. Beyond the horizon scale 

In the previous subsections, we have integrated the RG equations from the horizon scale (^=0) toward the smaller 
scale. It is also possible to solve them toward the larger scale with the same boundary conditions. The structure 
beyond the horizon scale is not observable at present but may give some relation with an inflationary scenario in the 
early Universe. In Fig. 7, we plot the scale dependence of (p), H and those fluctuations Ap and Ag in shear-free case 
with initial fluctuation Aq = 10~^'^. In the finite scale 10~^ < L < 10^, the flow stays in the vicinity of fixed point 
E and the power law behaviors of Ap and A^ toward larger scale are determined by Amax at fixed point E, while in 
larger scale L > 10^, the fiow converges to fixed point Q after H becomes negative and their power law behaviors are 
characterized by the eigenvalues of Wl^l at the fixed point Q. These eigenvalues are the same ones naively determined 
by the original scaling. This is consistent with the fact that the expansion ceases at the fixed point Q and the RG 
transformation reduces to a pure scale transformation. This property makes Ap constant around fixed point Q because 
of the cancellation between numerator and denominator in the definition of (|5^) (Fig. 7 (c)), while Ag diverges at 
the scale where H vanishes and turns negative from the definition ( ^ ) (Fig. 7 (d)). The declination of the slope in 
the log(p)-logL diagram (Fig. 7 (a)) is —2s which is the minimum eigenvalue of Wl^l associated with the fixed point 
Q. If we assign a fractal dimension D for this Universe, 13 = 3 — 2s because (p) reduces as e^^"''-'^ from a fractal 
structure while (p) oc e~^*^ by our scaling law. If we assign the value s = 0.5, the corresponding fractal dimension 
becomes D = 2. Because our model is a Newtonian cosmology, further analysis beyond the horizon scale would be 
irrelevant to the real Universe. 



IV. CONCLUDING REMARKS 

The renormalization group method has been applied to Newtonian cosmology in this paper and the scale dependence 
of the averaged observables on a fixed time slice including effects of fluctuations has been studied. The scaling 
assumption we proposed for the averaged observables is fulfilled in the case of either homogeneous distribution or 
non-analytic distribution such as a fractal, if we assume such a scaling at any position and at fixed time. Actual 
observation of the two point correlation function of galaxies may favor such a fractality at least in the finite scale 
range below 10 Mpc. Multi-fractality could be applied to explain the transition from finite fractal range at smaller 
scale to homogeneous distribution at larger scale . Our scaling analysis of averaged observables here might give a 
hint to the behavior of observables in such a fractal matter distribution, if it exists. 

To close our dynamical system and apply the RG method, here we have adopted an approximation of tidal force 
(LTA) and considered only fluctuations up to the second order. The higher order cumulants of averaged observables 
are ignored. We have found three robust flxed points of the RG equations. 

(1) The fixed point Q ("Quiescent" Universe) is stable in any directions toward larger scale. 

(2) The fixed point E ( "Einstein-de Sitter" Universe) is a saddle point. 

(3) The fixed point M ("Milne" Universe) is unstable in any directions toward larger scale. 

These fixed points do not change their stability even if we extend the "phase" space including fiuctuations as well 
as shear and tidal force. Any other fixed points are saddle, and therefore the scaling property of the global RG flow 
toward larger (smaller) scale is determined by fixed point Q (M). We find that the Universe asymptotically approaches 
fixed point Q or diverges toward infinitely regardless of the detail at the smaller scale. 

In order to find out the flow line which represents our real Universe, we have imposed the boundary condition flxed 
point E with tiny fiuctuations at the horizon scale. The inflationary scenario in the early Universe would imply such 
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a boundary condition, i.e. the inflation predicts the Einstein-de Sitter Universe at the horizon scale and beyond. We 
have solved the RG equation toward smaller scale from this boundary. Then we find that, toward smaller scale, the 
expansion H defined by (6')/3 monotonically increases from 2/3 to 1, and the density parameter = 8ttG{p) /3H'^ 
monotonically reduces from 1 to 0, provided the fluctuations of shear and tidal force are much smaller than the density 
fluctuations, which would be supported by the inflationary scenario |30| and observations of CMB. 

We have two important free parameters, s and Aq. These values might be fixed by some of the following arguments: 

(i) The scale dependence of the Mass-Luminosity ratio of astrophysical objects. 

(ii) A systematic decrease of the Hubble parameter H, if any, toward the larger scale. 

(iii) The scale dependence of the amplitude of the two point correlation functions. 

(iv) Some of the infiationary models predict a fractal-like structure beyond the horizon scale |3^ . 
Several comments related with our future work are in order. 

(a) We definitely need the general relativistic generalization of our work. This is because the large scale structures at 
high rcdshift and beyond the horizon scale would be adequately described only by general relativity. 

(b) In our analysis, we assume a scaling property of the averaged observables and their fluctuations. However we do 
not know whether such a solution is dynamically stable or not. To answer for this, we have to consider the dynamical 
stability of our system. 

(c) So far we have considered fluctuations of variables up to the second order. Though this would be sufficient when 
we consider the vicinity of the horizon scale, we ultimately need to incorporate higher order fluctuations as well for 
the consistency of the theory. 

(d) We have used the Local Tidal Approximation (LTA) for obtaining equations of motions for tidal force. By this 
method, we could avoid directly solving Poisson equation. This approximation would be justified upto mildly non- 
linear regime. However it may not be valid in the fully non-linear regime, especially in smaller scales, where the 
fluctuations grow Ap 3> 1. In order to treat the smaller scales properly, we need to improve LTA in our future study. 
We would like to report the analysis on these issues in the near future. 
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Table 1. The list of three fixed points and their eigenvalues in the case with neither shear and tidal force nor 
fluctuations. 
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Table 2. The list of 7 fixed points and their eigenvalues in the case without fluctuations in LTA. The seventh fixed 
point [7] is not an isolated point but a closed loop parameterized by one parameter 77 as cr^ = ^ cos 77, crl = | sinr], 

= |(cos2?7 — 3 COST/), i?* = — ^(sin2r7 + 3sinr7), where < 77 < 27r. Ai, A2, A3 are three real solutions 
(-1 < Ai < A2 < 0,A3 > 1) of A^ - A-4(l-cos3r?)/27 = 0. 
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Table 3. The list of fixed points and their eigenvalues in the case without shear and tidal force including 
fluctuations. New fixed points (8) ^ (10) are unphysical because (p^)* and/or (^^^)* are negative. 
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Table 4. Three important fixed points and their eigenvalues in the case with shear and tidal force in LTA as well as 
fluctuations. 
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APPENDIX A: THE BASIC EQUATIONS FOR THE CASE WITH FLUCTUATIONS 



To be complete, here, we present a full set of our basic equations, which have been used in our renormalization 
analysis for Newtonian cosmology. Those include fluctuations up to second order cumulants as well as shear and tidal 
force. 



(A) The RG equations for averaged vEiriables 

d{p) 



di 

m 

d£ 

d{<J+) 
di 

d{<T-) 
d£ 

d{E+) 

de 

d{E-) 
d£ 



= S 



= s 



= s 



= s 



lie) - l] - {a+r + {a_r + - ^{0a+), - (4), + (al). 



m - 2) {E+) + AnG{p){a+) + AnG{pa+), 
- 2) {E_) + 4nG{p){a-) + 4nG{pa-), 



(B) The RG equations for second order cumulants 

d{p'}o 



de 

d{6^)c 

de 

d{p9), 

de 

d{p(T+)c 

de 

d{pa-)c 

de 

d{pE+), 

de 

d{pE_)^ 

de 



= 2S 



= 2S 



(W-2)(p^)e + (p)(p^)c 



1 {9^), + 47rG(p0)e + 12 {{a+){9a+), + {a_){9a_),) 



= S 



= S 



= s 



-{9} - 3j {p9), + 47rG{p% + {p){9\ + 12 {{a+}{pa+), + (a_)(pa_)e) 
^{9) - 2{a+) - 3^ + 2{a.){pa.), + |((7+)(p0)e + {p){9a+), + {pE+), 

'^{9) + 2{a+) - a") (pa_)c + 2(a_)(pa+)e + \{a-){p9), + {p){9a.), + {pE_), 



2 {{9} - 2) {pE+), + {E+)(p9), + {p){9E+), + AnG {{a+){p^), + {p){pa+),) 
2 {{9} - 2) ipE-), + {E.){p9)c + {p){eE-), + 47rG ((a_)(p^)c + {p){pcT-)c) 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 



(A7) 

(A8) 

(A9) 
(AlO) 

(All) 

(A12) 

(A13) 
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d£ 



d{0E_), 
d£ 



d{a%)c 
di 

d{al)c 
dl 

d{a+a-)c 
di 



d{El), 
di 

d{E'i), 
di 

d{E+E_), 
di 



d{a+E+)^ 
di 



S 



2 {^{6)- {a+) - l) {ea+)c + 2{a-){6a-), + \{cj+){6% + {6E+), + 47rG(pa+)c 



+ 12 ((a+)(a2), + (a_)(a+a_)e) 



: S 2 {^{e) + - 1^ {6a_), + 2(a_)(^a+)e + \{cJ-){e^)o + {0E_), + 47rG(pa_), 
+12 ((a_)(a2), + (a+)(a+(7_)c) 



^{9) - 3) {eE+), + {E+){e\ + 47rG {{p){ea+), + + {pE+),) 



+12 {{a+){cj+E+), + (a_)(a_i;+)c) 



^-{e)-i] {eE_), + (£;_)(02^c + 47rG ((p)(0a_)c + + {pE_),) 



+ 12 ((a+)(<7+i;_)e + ((T_)(<7_i;_)e) 



= s 



2 - 2(a+) - 1^ {al), + ^(a+)(0a+)c + 4(a_)(a+a_)e + 2(a+i?+)e 

2 0(0) + 2(a+) - 1^ {al), + ^(a_)((?a_)e + 4(a_)((7+a_)e + 2{a.E_), 
2 0(0) - 1) (a+a_)c + I ((a+)(0a_)c + (a_)(0a+)c) + 2(a_) ((4), + (a^),) 

+(a+i;_)c + (a_i;+), 

2S {{9) - 2){El), + {E+){eE+), + 47rG((p)(a+i;+)c + (a+)(pi?+)c) 
{{9) - 2){El), + {E_){eE_), + 47rG((p)(a_£;_)c + (a_)(p£;_)c) 
2 {{6) - 2) {E+E_), + {E_){eE+), + {E+){9E_), + 47rG(p) ((a+i?_)c + (a_i?+)c) 



= 25 



= 5 



+47rG {{a+){pE_), + (a_)(p£+)c) 



= 5 



0(0) - 2(a+) - 3^ {a+E+), + 2(a_)(a_£;+)e + (£;+)(0a+)e + {El), 



+-{a+){eE+), + AttG {{p){al), + {a+){pa+). 
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dl 



S 



{6) - 2{a+) - 3 ((T+£;_)c + 2(<T_)(a_£;_)c + {E-){ea+)c + {E+E. 



+-{a+){eE_), + 47rG((p)(a+a_)e + (a_)(p(7+)c) 



d{a-E+), 
dt 



= S 



-{9) + 2{a+) - 3 ) (a_^+)e + 2{a-){a+E+), + {E+){ea-), + {E+E. 



+-{a-){eE+), + A7TG{{p){a+a.), + {a+){pa-U) 



d{a-E_)c 



= S 



(0) + 2{a+) - 3 (a_£_)c + 2{a-){a+E_), + {E_){ea-), + {El), 



+-{a.){eE_), + 47rG {{p){al). + (c7_)(pa_)c) 
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FIG. 1. The RG flow toward larger scale on the H-Q plane in the case with neither shear nor fluctuations. Here we define 
the Hubble parameter H and the density parameter Q as H = {0)/3 and CI = 87rG(p)/3i/^, respectively. The three important 

fixed points exist: (1) fixed point Q (Quiescent Universe), (2) fixed point E ( "Einstein-de Sitter" Universe) (3) fixed point M 
("Milne" Universe). The RG fiow converges to the fixed point Q or escapes to infinity. The fixed point seems to be a line in 
the H-Cl plane, but it is a point that all physical variables vanish in the "phase" space of {F) = ((p), {0}). fl turns out to be 
some finite value, which depends on the limit of (p) —^ and {6} —^ 



FIG. 2. The scale dependence of (a) logfJ, (b) H, (c) Ap and (d) Ag in shcar-frcc case with s = 0.5 and 
Ao = (1) lO^S (2) 10"^ -\ (3) 10"^ (4) 10"^ ^ We solve the RG flows from the fixed point E toward smaller scale. 
L = exp£ represents the ratio of the scale of averaged region to the horizon scale. Every flow monotonically converges to M as 
H becomes laxger {H /") and Q, gets smaller {Q, \) toward smaller scale. The critical scale where Q. and H deviate from the 
EdS values depends on the magnitude of Aq. 



FIG. 3. The s dependence of the RG flows of (a) log Q., (b) H , (c) Ap and (d) Ag from the fixed point E at horizon scale in 
shear-free case. Here we depict the cases of Ao — 10^^ and s = (1) 0.3, (2) 0.5, (3) 0.7, (4) 0.9. Q, begins to decline beyond 
a critical scale, a few decades smaller than horizon scale, which depends on the magnitude of Aq. The slope at smaller scale is 
determined by aJ^Ij^ = s/(l — s) at the fixed point M. The slope for the case (2) is equal to 1. H also changes from the value at 
E to the value at M beyond the above critical scale. As the value of s increases, the critical scale shifts toward smaller scale. 



FIG. 4. The RG flow toward smaller scale on the H-Q, plane in the case that s = 0.5 and Acr^.o = A_e^,o = 10 ^ . As for the 
density perturbations, we choose Ao = (1) 10"° '^®, (2) 10"\ (3) lO"^'^^, (4) lO"^'^. If both A<^^,o/Ao and Ai5_^,o/Ao are 
less than Rcr ^ 0.1495 ((1), (2)), every flow converges to M as increases {H y) to 1 and Q decreases (fi \) to toward 
smaller scales. On the other hand, if Aa^,o or Ae^,o are larger than the above critical value Rcr ((3), (4)), the flow heading 
for M comes to turn around and go away to inflnity. 



FIG. 5. The RG flow of (a) logO, (b) H, (c) \og{{a±)/H) and {d)\og{{E±) / H'^) from the flxed point E toward smaller scale 
in the case that s = 0.5, A<,^,o = Ab_^,o = 10"^ and Ao = (1) 10-° '"^, (2) lO'S (3) IQ-^ '^^ , (4) IQ-^-^. In the case of (1) 
or (2) (Act^,o/A(). A_E^.o/Ao<-Rcr ~ 0.1495), the RG flow converges to M toward smaller scales as the case without shear and 
tidal force (sec Fig. 2). In the case of (3) or (4) (A^^.o/Aq, A_B^,o/Ao>_Rcr)! however, the RG flow does not converge to M 
but rather diverges to infinity. We also depicted the flow of fluctuations; (c) log Ap, (f) logAg, (g) log A(t+ and (h) logA_B+ 
for the same initial values as the above. In the case of (1) or (2), all fluctuations but Ap damp off as the RG flow converges to 
M, while Ap approaches the finite value. In the case of (3) or (4), all of them diverges to infinity. 



FIG. 6. The s dependence of the RG flows from fixed point E of (a) logfi, (b) H, (c) log({cr±)/H) and {d)log{{E±) / H'^). 
Here we depict the cases of s = (1) 0.3, (2) 0.5, (3) 0.7, (4) 0.9 with the boundary conditions, Ao = lO"^'^ and 
'^CT+.o ~ '^-B+.o ~ 10^^ "' at horizon scale. The results for log SI and H arc quite similar to Fig. 3, if the flow converges 
to the fixed point M. The slopes of (cr±) and {E±) at larger and smaller scale are determined by at fixed points E and 
M, i.e. — 4s/(3 — 2,s) and 2s/(l — s), respectively. As s increases, the turning point of the slope shifts toward smaller scale and 
the slope becomes less steep. We also depicted the flow of fluctuations; (e) logAp, (f) logAg, (g) logA^^ and (h) logAf;, 
for the same range of s. Ao = 10~^'® and Act_,.,o = ^b+,o = 10~^'® at horizon scale in the case of LTA. We find that Ap 
= constant near the flxed point M, while it decreases linearly near the fixed point E, which behavior is determined by the 
eigenvalue —4s/ (3 — 2s). 

FIG. 7. The scale dependence of (a) log(p), (b) H, (c) logAg and (d) logAp beyond horizon scale for s = 0.5 with the 
boundary conditions, Aq = 10~^'® at horizon scale in shear-free case. In the flow departing from the fixed point E, (p) 
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approaches zero toward both smaller and larger scales, for which the slopes at larger and smaller scale are determined by X^\^ 
at fixed points, Q and M, i.e. —2s and 2s/(l — 2s), respectively. H also vanishes at larger scale, after it becomes negative 
around L = 10®, which causes the divergence of Ag around the scale. The slope of both Ap and Ag around the flxed point E 
are determined by a|^!jj = — 4s/(3 — 2s) toward smaller scale and AmLx = 6s/(3 — 2s) toward larger scale at the point. 
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